Electrooptics of graphene: field-modulated reflection and birefringence 
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The elecrooptical response of graphene due to heating and drift of carriers is studied theoreti- 
cally. Real and imaginary parts of the dynamic conductivity tensor are calculated for the case of 
effective momentum relaxation, when anisotropic contributions are small. We use the quasiequilib- 
rium distribution of electrons and holes, characterized by the effective temperature of carriers and 
by concentrations, which are controlled by gate voltage and in-plane electric field. The geometry 
of normal propagation of probe radiation is considered, spectral and field dependences of the re- 
flection coefficient and the relative absorption are analyzed. The ellipticity degree of the reflected 
and transmitted radiation due to small birefringence of graphene sheet with current have also been 
determined. 

PACS numbers: 78.20.Jq, 78.67.Wj, 81.05.ue 



I. INTRODUCTION 

Study of electrooptical response both of bulk semicon- 
ductors and of heterostructures (see Refs. in [l[ and 0, 
respectively) is a convenient method to characterize these 
materials. Such a response is used to modulate both 
the intensity of radiation and its polarization. As it was 
demonstrated more than 30 years ago [1,0] (see also Sect. 
17 in [5|]), the main contribution to the elecrooptical re- 
sponse of narrow gap semiconductors is caused by the 
modulation of the interband transitions, both virtual and 
real one, under heating and drift of nonequilibrium carri- 
ers. The electrooptical properties of two-dimensional car- 
riers in heterostructures have also been studied [6( . Since 
graphene is a gapless semiconductor with linear energy 
spectrum 0] , the direct interband transitions in graphene 
are allowed with the characteristic interband velocity 
vw = 10 s cm/s, which corresponds the Weyl-Wallace 
model [H, degenerated over spin and valleys. Therefore, 
optical properties of graphene should be modulated es- 
sentially by temperature and carriers concentration Q 
and these dependences were studied recently, [lfj The 
applied electric field not only changes carriers tempera- 
ture and concentration, but also causes the anisotropy of 
distribution due to carriers drift [ll|, [T§]. Therefore, the 
birefringence effect can be essential for radiation prop- 
agating across a graphene sheet with current. To the 
best of our knowledge, no measurement of electrooptical 
response was performed until recently, and a theoretical 
study of these phenomena is timely now. 

The results obtained below are based on the tensor of 
dynamic conductivity, determined by interband transi- 
tions of non-equilibrium carriers. This tensor is deter- 
mined by Kubo formula in collisionless approximation 
(3) HI with the use of weakly aniso trop ic distributions of 
electrons and holes calculated in [lj, [lj]. The case of 
normal propagation of the incident (in), reflected (r), 



and transmitted (t) waves of probe radiation (see Fig. 1) 
is studied, and the reflection and transition coefficients, 
controlled by carriers heating conditions, are obtained. 
It is demonstrated, that the heating level dependence on 
applied field, temperature of phonons, and sheet charge, 
controlled by gate voltage V g , can be verified from elec- 
trooptical measurements. Moreover, graphene is to be 
considered due to carriers drift as an uniaxial plane, and 
the elliptically polarized r- and i-waves appear under 
linear polarization of m-radiation, if 9 ^ or 7r/2, see 
Fig.l. Due to an effective relaxation of carriers momenta 
the distribution anisotropy and the induced birefringence 
are small, but a high sensitivity of polarization measure- 
ments enables one to determine drift characteristics of 
nonequilibrium carriers using a field- induced Kerr effect. 




'Electronic address: ftvasko@yahoo.com 



FIG. 1: Schematic view on incident (in), reflected (r), and 
transmitted (t) radiation for the case of normal propagation 
through the graphene sheet with applied electric field E. An- 
gle 9 defines the polarization direction of in- wave while r-, 
and t-contributions are elliptically polarized. 



The consideration below is organised in the following 
way. In Sec. II we present both the basic equations for the 
complex tensor of dynamic conductivity, and the electro- 
dynamics equations for the uniaxial graphene sheet on a 
substrate. Numerical results, describing the electromod- 
ulation spectra and Kerr effect, are discussed in Sec. III. 
The concluding remarks and the list of assumptions are 
presented in the last Section. In Appendix, the dynamic 
conductivity of the undoped graphene is considered. 
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II. BASIC EQUATIONS 

The description of graphene response on the probe in- 
plane electric field E w exp(— iu>t) is based on the con- 
sideration of the high-frequency dynamic conductivity 
and on the examination of the electrodynamics problem 
for propagation of such a field through graphene sheet. 
When performing these calculations, we take into account 
a modification of interband transitions due to carriers 
heating, and an anisotropy of response due to carriers 
drift. 



A. Anisotropic dynamic conductivity 

Contribution of the interband transitions of non- 
equilibrioum carriers with the distribution function /; p , 
into the response at frequency u> is described by the dy- 
namic conductivity tensor a a p{uj) given by Kubo for- 
mula: 



A(ev w ) 2 



dip - ffp) 
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Here |Zp) and Ei p are the state and energy of Ith band 
(c- or v-) with the 2D momentum p, and A — > +0. We 
also use the velocity operator vw<? and the normalizing 
area L 2 . The expression (1) is written in a collision- 
less approximation lot >> 1 (r is the relaxation time), 
when intraband transitions are inefficient. In this case 
the density matrix, averaged over scattering, should be 
used in Kubo formula, and a a p(uj) appears to be written 
through the stationary distribution function /; p . [5[ Due 
to effective momentum relaxation the anisotropy of non- 
equilibrium electrons and holes distributions is weak and 
the expansion 



fip — fip 



A// p 2) cos2v? 



(2) 



should be substituted into Eq. (1). Here ip angle deter- 
mines orientation of p and Afj^ oc E r , where E||OX is 
a dc field applied. The linear in E contribution can be 
omitted from a a p{ui) in the case, when the small spatial 
dispersion, responsible for the radiation drag by current 
(see Ref. 14), is neglected. Thus, with an accuracy of 
the contributions of oc E 2 order, tensor (1) is determined 
by the transition matrix elements: 



i<ipi^.»i-ip>r = 1/2, 

-1/4, 
1/4, 



cos 2ip | (lp \a x \ - lp) | 



(3) 



cos2<£ |(lp \a y \ - lp)|" 
where overline means the averaging over angle. 



Since the non-diagonal components of tensor (1) van- 
ish, the XX- and YY -components of dynamic conduc- 
tivity: 



Act,, 



(4) 



describe the response of graphene sheet with the field- 
induced uniaxial anisotropy. Further, we substitute Eqs. 
(2) and (3) into (1), and we use the electron-hole rep- 
resentation, when f cp — > / ep , and f vp — > 1 - fh,- P , 
see [ll|. It is convenient to separate the contributions 
of the undoped graphene and of the free carriers (elec- 
trons and holes) into the isotropic part of conductivity, 
~g u and , so that er w = er w + crjy . The first contri- 
bution is discussed in the Appendix. Separating the real 
and imaginary parts of , and using the energy conser- 
vation law, one obtains the frequency-independent result 



Rcer 



/47i. The imaginary contribution into a u is 



given by the phenomenological expression (A. 2), which 
is written through the fitting parameters corresponding 
to the recent measurements. |15| 

The electron- hole contributions to the isotropic part, 

(c) 

<Tu , and the anisotropic addition, Au w , are written as 
follows: 
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The real parts of conductivity given by Eqs. (5) are 
expressed directly through isotropic distribution (2) at 
the characteristic momentum for interband transitions, 
hu;/2vw, according to: 
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The imaginary parts of affl, and Act^, given by Eq. (6) 
are transformed into: 



Im 
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and the principal value integrals here should be calcu- 
lated numerically. 
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Below, we restrict ourselves to the case of quasielastic 
distribution of non-equilibrium electrons and holes (k = 
e, h) with effective temperature T c and chemical potential 



f kp = {exp[(v w p - ptk)/T c ] + 1}~ 



(9) 



The dependences of distribution (9) on electric field E, 
temperature T, and gate voltage V g are presented in [ill . 

Il3j ]. For the anisotropic addition Aff^} in the case of 
momentum relaxation through elastic scattering we use: 



2^ 2) dp 



(i) 

\_PVp 



dfkp 
dp 



(10) 



For the case of short-range scattering on static defects 

(12) 

the relaxation rates v p ' are proportional to the density 

of states, so that Vp — Vdp/Ti + vq, and Vp = 2vp + 
Uq. Here Vd is a characteristic velocity, that determines 
an efficiency of momentum scattering, [l6[ and vq is a 
residual rate, which describes the scattering process for 
low-energy carriers. 



B. Electrodynamics 

For normal propagation of probe radiation, the Fourier 
component of the field, E W2 , is governed by the wave 
equation: 



dz 2 



Attuj . 



0, 



(11) 



where e z is dielectric permittivity. In this paper we ex- 
amine the case of graphene on the thick SiC>2 substrate, 
when e z <o = 1, and e z> o = e. The induced current 
density in (11) is localized around z — plane, so that 

+o 

jwz ~ iTuzEu^o, while / dza wz = <j u is determined 

-o 

through the dynamic conductivity tensor, being exam- 
ined above. Outside the graphene sheet the solution of 
(11) can be written as: 
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Here the wave vectors fc w = lu/c and k u — ^/euj/c, and 
the amplitudes of incident, reflected, and transmitted 
waves (E( m \ E' r ), and E^ respectively) were intro- 
duced. This amplitudes are governed by the boundary 
condition, 



dE„ 



dz 



, 47T 

+ ikuj — cr w E W2= o = 0, 
c 



(13) 



which is obtained after integration of Eq.(ll) over z 
through the graphene layer (—0 < z < +0). The sec- 
ond boundary condition is the requirement of continuity: 



E, 



E, 



Taking into account the diagonality of cr^ tensor, we 
get the solutions from the boundary conditions as follows: 



E (t) 



l + A a (u) 



1 + A a (uj) 



(14) 

where factor A a (w) = \[i + Ait<j aa {uj) / c was introduced. 
After substitution of Eqs. (12) and (14) into the general 
expression for Poynting vector S = c 2 Re [E x rotE*] /8ir, 
we obtain the incident, reflected, and transmitted fluxes 
Sim S r , and S t respectively, which are parallel to OZ. 
After multiplication of Eq. (13) by Ejf, we get the rela- 
tion between these fluxes as follows: 

S in = S r +S t +-Re{W t -a UJ --E lt ), (15) 

where the last term describes absorption. 

The polarization characteristics of r-, and i-waves are 
determined by solutions (14). It is convenient to present 
them in complex form E a = £ a exp(z-0Q,), where E a , and 
ifj a give the amplitude and phase of a-component of the 
field respectively. At 8 = 0, or at 8 = 7r/2, when the 
response is described by a xx , or by <7 yyi the linearly 
modulated r-, and i-waves occur. For other 8, the re- 
flected and transmitted waves are elliptically polarized. 
The ellipticity degree e(oj) is determined by the phases 
difference between X- and y-components of the field, 
Atjj = tp x — tjjy, see the general expressions in Ref. 17. 
Under weak anisotropy, with the accuracy up to first or- 
der in Acr^, we get s(oj) — A%p/2. 



III. RESULTS 

Now we examine the spectral and polarization char- 
acteristics of the electrooptical response. We study the 
reflection, transmission, and relative absorption coeffi- 
cients, determined as R^g = S r /Si n , T u g = St jSi n , 
and £ u = Re (Ej • (7 W • E t ) /25 m , respectively, [3 as 
well as the ellipticity degree, e(ui), for the case of weak 
anisotropy. The final expressions for the coefficients un- 
der consideration are obtained with the use of complex 

conductivities ajj \ and Aa^, given by Eqs. (7), and (8), 
and they depend both on Hu>/T c , and on carriers concen- 
tration. In Fig. 2. we plot these dependences and one 
can see that the response modify essentially with tem- 
perature and concentration. The smallness of anisotropic 
additions is determined by dimensionless factor 



F = 



zEhv 



w 



2T 2 



Vd 



(16) 



uiz=+0- 



which arises from oc E 2 contribution to the distribution 
function (10). Note also, that ImAer^ depends weakly 
on the cutting parameter {hv^vw) / {T c Vd), taken below 
as 0.1. 
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FIG. 2: Spectral dependences of <r£, c ' (a) and Acr^ (b) for 
intrinsic graphene with / P =o =0.5, 0.3, and 0.1. Solid and 
dashed curves correspond real and imaginary parts of con- 
ductivity, respectively. 



A. Reflection and absorption 



0.75- 



0.50 



0.076 



0.075 



£ = 30 V/cm 
20 V/cm 
-10 V/cm 
OV/cm (a) 



0.05 0.10 0.15 0.20 





£=30 V/cm (c) 
.20 V/cm 
10 V/cm 



0.05 



0.10 0.15 

Ao(eV) 



0.20 



FIG. 3: Spectral dependencies of relative absorption (a), 
reflection (b) and differential reflectivity (c) for intrinsic 
graphene at 77 K and at different electric fields, E (marked). 
Solid and dashed curves in (b) are plotted at E =0 and 30 
V/cm, respectively, for e m =60 meV (1), 80 meV (2), and 100 
meV (3). 



For the examination of R^e, and £ w g it is convenient 
to separate the isotropic and ^-dependent contributions, 
so that 



R U) +AR u cos2e, ( u 



■A&,cos20, (17) 



where the small (of the order of Acr^/er^) anisotropic 
additions, proportional to cos 29, have been separated, 
see Fig. 1. The coefficients in Eq. (17) are written below 
through <j w , Acr w and the factor A u = yfe + Ama^jc. For 
the isotropic parts of reflection, and relative absorption 
coefficiens we get [1]: 



Ru. 



1- A 



1 + A U 



IQn Recr w 



(18) 



so that these characteristics depend on T, E, and V g . 

Spectral dependences of the relative absorption, 
(,ui/£,max, the reflection, R^, and the differential reflectiv- 
ity (5R/R) U = (Ru - R^ q) )/Ri eq) , are plotted in Fig.3 
for intrinsic graphene at T = 77 K and different electric 
fields (the data for T c and carriers concentration were 



used from Ref. 11). Here £, max is the maximum value 
of relative absorption for high frequences, when the free 
carriers contribution is unessential. One can see, that 
due to the increase of average energy of carriers with the 
increase of E the absorption increases at high frequencies 
and decreases for the low ones. The relative change of £ w 
is reasonably large, and for Tilu ~ T c it can be measured 
directly. At the same time, the reflection coefficient de- 
pends on field in more weak way, see Figs. 3b, c, and 
(8R/R) U can be 10 -2 in THz spectral region; in near-IR 
spectral region it decreases down to value < 10~ 3 . Note, 
that for hui <0.1 eV R u increases essentially (at high fre- 
quencies R^ ^0.075) due to the contribution of the first 
summand of Eq. (A2). Fig. 3b presents the dependence 
of R u on phcnomcnological parameter e m ; £ w depends 
weakly on this parameter. 

The dependences of R^, and (5R/R) U on doping 
level are presented in Fig. 4. The data for the room tem- 
perature are presented for V g =3 V and 10 V, which 
correspond the difference between electron and hole con- 
centrations 1.65 x 10 11 cm -2 and 5.5 x 10 11 cm -2 , respec- 
tively. Similarly to field dependences at T =77 K (see 
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FIG. 4: The same as in Fig. 3 for doped graphene at room 
temperature and different V g (marked). Solid and dashed 
curves correspond E =0 and 30 V/cm, respectively. 



Fig. 3), with the increase of V g (the doping level oc V g ) 
the response moves towards the high energy region. The 
dependences on the level of heating (the applied field E), 
and on carriers concentration (the gate voltage V g ) corre- 
spond the measurements of spectra for different tempera- 
tures and V g , see [l(| • For the range of parameters under 
examination the field modulation of £ w is of 20-^50 % 
order up to mid-IR spectral region. These modifications 
should be observed rather easily. The carriers contribu- 
tion into reflection increases as well: at hu> >0.1 eV the 
decrease of R u occurs, which almost does not depend on 
e m ; in this case the shape of the differential reflectivity 
is similar to low temperature case, with the shift into the 
high energy region. 

Later we shall examine the anisotropic contributions 
in Eq.(17), which are proportional to Aa u . Such a con- 
tribution into reflection coefficient is given by: 
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and the addition to relative absorption takes form: 

'27rAcr w /c 



(19) 



A£,= 



167r RetJcj 



1 + A 



+ c.c. 



0.9 



\ 0.6 
^ 0.3 
0.0 
^ 0.9 
51 0.6 
0.3 



0.0 



\ £ = 30 V/cm 




\\ / 20 V/cm 


(a) 


\ Y / 10 V/cm 


\\\/ / OV/cm 




X. £ = 30 V/cm 




\/ 20 V/cm 


(b) 


\ \/ 10 V/cm 


\\\ xS<^ V/cm 







10 20 30 40 

hco(meV) 

FIG. 5: Spectral dependences of anisotropic contributions 
to relative absorption, A^/F (a), and to reflectivity, 
ARu,/(RujF), (b) for intrinsic graphene at T =300 K and 
at different electric fields, E. 
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Spectral dependences for anisotropic contributions to 
the relative absorption and reflectivity, A^/F and 
AR LJ /(R U1 F), are shown in Figs. 5 (a) and (b). One 
can see, that in the range of fields under examination the 
parameter F given by Eq. (16) does not exceed 0.05, so 
that A£ w and AR^/R^ are of 10~ 4 order for the mid-IR 
spectral region (~ 0.1-^0.2 eV) and the response increases 
up to ~ 10~ 2 in THz spectral region. The anisotropy of 
such order of value can be analyzed with the use of the 
modulation methods. 



B. Kerr effect 

Besides the cases of parallel or transverse orientation of 
the probe radiation polarization with respect to the drift 
direction (i. e. at 6 ^ 0, tt/2), the reflected and trans- 
mitted fields are elliptically polarized. The maximal Kerr 
effect occurs if the z-wave is polarized along 9 = ir/4, and 
below we consider this case only. In the approximation 
of the weakly anisotropic distribution (2) the ellipse ori- 
entation does not differ essentially from 6 ~ ir/4, and the 
ellipticity degree can be written as (l~7| : 



e(u) = sin/3Re$(w) — cos/3Im$(w). 



(21) 



Here the (3 angle, and the complex function <S>(w) can 
be expressed through the difference of the phases of r-, 
and /-waves (see the end of Sec. II). The smallness of the 
ellipticity is determined by the relation $(w) oc F, while 
the (3 angle is not small. 



6 



fen 



0.08 



0.04 



0.00 





-0.04 




0.08 : 




0.06 : 




0.04 : 


to" 






0.02 : 



£ = 30 V/cm 

20 V/cm 
10 V/cm 
V/cm 





£ = 30 V/cm 
20 V/cm 
10 V/cm 
V/cm 



25 50 75 100 

/z£)(meV) 



125 



FIG. 6: Spectral dependences of ellipticity degrees of r- and t- 
waves [panels (a) and (b), respectively] for intrinsic graphene 
at T =300 K and different electric fields, E. 



For the reflected wave the $(w) function is given by 
the expression 



$ r (w) = 



47rAcr„ 



1 + A 



1 - A 



c(i + A 

while the /3 angle is introduced through the relation: 

21mA 
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l-K 
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Similarly, for the transmitted wave, (20) is expressed 
through the function: 



2wAa, 



11 + A 



c(l + A) 2 ' 
and the [3 angle is given by the expression: 

ImA 

tan [3t 



1 + Re A 
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Substitution of these expressions into Eq. (21) gives the 
ellipticity degrees for r- and t- waves, £ r (u>) and et(to). 

Spectral and field dependences of e r (uj) and e t {<-o) are 
shown in Figs. 6a and 6b for intrinsic graphene at T =300 
K. In mid-IR spectral region e rt t(uS) decreases with ui and 
for Ku> ~ 0.1-f-0.2 eV the value of ellipticity degree does 
not exceed ~ 10~ 4 at F ^0.05. In THz spectral region 
e r ,t{w) increases up to ~ 5- 10 -3 , wherein the direction of 
rotation for the reflected wave changes at hu> ~25 meV. 
Such value of ellipticity degree can be detected by mod- 
ulation methods only. However, in stronger fields, when 
the distribution function is strongly anisotropic, [121 ] the 
ellipticity degree can increase essentially. 



IV. CONCLUSIONS 

Summarizing the consideration performed, we have ex- 
amined the graphene electooptical response due to the 
interband electron transitions under the carriers heat- 
ing and drift. It was found, that an essential modula- 
tion of the reflection, and the relative absorption take 
place starting from the field strength ^30 V/c at liquid 
nitrogen and room temperatures (with the increase of 
field, the modulation should increase essentially). Due 
to current-induced birefringence of graphene sheet the 
weak ellipticity of the reflected and transmitted radia- 
tion arise. 

Next, we list and discuss the assumptions used in our 
calculations. First, the dynamic conductivity tensor (1) 
is written in collisionlcss approximation. For the case 
of short-range scattering, when uj >> v^p^/h, one ar- 
rives to the condition Vd/2vw <C 1 and the collisionlcss 
approximation is not valid for a strongly disordered ma- 
terial. Also, the interband response of a pure graphene 
is described with the use of the phenomenological ex- 
pression (A. 2) and a low-frequency restriction for this 
approximation is not clear. 

Second, the quasiequilibrium distribution of carriers 
(9) was used for the numerical estimation of electroop- 
tical response. This means the assumption of an effec- 
tive intercarrier scattering. The complete description of 
the carriers heatin g un der such conditions had not been 
performed yet [TlTTl9j. However, the approximation (9) 
gives a good estimation for the response magnitude, and 
the peculiarities of spectral dependences enable us to 
determine the contributions of the different relaxation 
mechanisms. Similarly, the use of short-range scattering 
model in the drift-induced contribution (10) gives the es- 
timation for optical anisotropy magnitude, and the spec- 
tral dependences peculiarities contain information about 
the momentum relaxation mechanism (despite the short- 
range scattering can be treated as a dominant one within 
the phenomenological description of momentum relax- 
ation, [Tfjj the microscopic mechanism have not been ver- 
ified until now (20j). 

Third, we have examined the heating of carriers with 
low energies only, (the results for E <30 V/cm have been 
presented) , when the essential electrooptical response oc- 
curs in THz spectral region only. With the increase of 
field (up to tens kV/cm, see [12|]) the electrooptic effect 
increases and shifts into near-IR spectral region. The the- 
oretical approach developed here can be applied for this 
region as well, however the calculation of the distribution 
of hot carriers for this case have not been performed yet. 

Forth, the case of graphene on a thick substrate have 
been examined. The consideration of the interference 
effects for graphene, placed on substrate of limited thick- 
ness, needs more complicated calculations (however, an 
accuracy of measurements can increase for near-IR spec- 
tral region [21]), and is beyond the frame of this paper. 
And the last, we have limited ourselves to the examina- 
tion of the geometry of normal propagation of radiation 
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only. The study of the response dependence on the angle 
of radiation falling gives additional experimental data, 
however it is more complicated and needs special treat- 
ment. 

Finally, the results obtained demonstrate, that the 
electrooptical response due to heating and drift of car- 
riers is large enough, and it can be measured. Because 
of strong dependence of the response on the applied 
field, temperature, and gate voltage, these measurements 
can give an information on relaxation and recombination 
mechanisms. In addition, the electrooptical response of 
graphene can be applied for modulation of intensity and 
polarization of radiation in THz and mid-IR spectral re- 
gions. 



where the real and imaginary parts of conductivity have 
been separated. The direct integration with the use of the 
energy conservation law gives the frequency-independent 
real part of (Al): Rea = e 2 /4h. The imaginary con- 
tribution into a^, appears to be divergent at p — ► oo; 
moreover Tma^ cx vwPm/hu, where p m is a cut-off mo- 
mentum. On the contrary to the case of bulk ma- 
terial, @ this cut-off appears to be too rough for the 
description of the response in graphene. It is convenient 
to approximate 7ma w , by separating the terms cx cj -1 , 
and cx u>, which correspond the contributions of the vir- 
tual interband transitions, and of ions, correspondingly. 
As a result, we get: 



APPENDIX: RESPONSE OF UNDOPED 
GRAPHENE 

The dynamic conductivity for the case of undoped 
graphene is described by Eqs. (1) and (3) after replace- 
ment of f vp by 1 and of f cp by 0. As a result we get the 
expression: 



2{ev w ) 
uiL 2 



E 



ir5(hto — 2vwp) 



iV 



Tloj — 2vwP 



(A.l) 



Imov 



h \huj 
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where the characteristic energies, s m , and £j, have been 
introduced. The comparison of the response, described 
by Hut with the recent measurements of the graphene 
optical spectra, yields: e m ~0.08 eV, and Ej ~ 6.75 eV. 
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